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We propose a QCD axion model where the origin of PQ symmetry and suppression of axion 
isocurvature perturbations are explained by introducing an extra dimension. Each extra quark- 
antiquark pair lives on branes separately to suppress PQ breaking operators. The size of the extra 
dimension changes after inflation due to an interaction between inflaton and a bulk scalar field, which 
implies that the PQ symmetry can be drastically broken during inflation to suppress undesirable 
axion isocurvature fluctuations. 


I. INTRODUCTION 

The null result of the electric dipole moment for neu¬ 
tron implies that the CP phase in the strong sector is 
smaller than of order 10” 10 [1]. An elegant mechanism 
to explain the smallness of the strong CP phase is the 
PQ mechanism [2, 3]. The strong CP phase is promoted 
to a dynamical field, called axion, which is associated 
with a spontaneously broken anomalous symmetry [4], 
The axion obtains a periodic potential through nonper- 
turbative effects [5, 6] and starts to oscillate around its 
minimum at the QCD phase transition. The strong CP 
phase is dynamically cancelled by the vacuum expecta¬ 
tion value (VEV) of the axion. In addition, the axion 
oscillation behaves like the cold dark matter (DM) and 
can explain the observed amount of DM [7-9]. However, 
the PQ symmetry is an anomalous global symmetry, so 
that it is not trivial that Lagrangian respects that sym¬ 
metry in a sufficiently high accuracy [10-13]. In fact, in 
the literature they pointed out that the PQ symmetry 
may be broken by quantum gravitational effects, so that 
severe fine-tunings may be required to solve the strong 
CP problem by the PQ mechanism [14-16]. 

In Refs. [17, 18], Izawa, Watari and one of the present 
authors proposed a model that has an accidental PQ sym¬ 
metry by introducing an extra dimension. 1 Each extra 
quark-antiquark pair lives on branes separately so that 
PQ breaking operators are suppressed exponentially by 
the size of the extra dimension. The extra quarks are 
charged under QCD and a bulk hypercolor that is as¬ 
sumed to be confined at an intermediate scale of order 
10 12 GeV. Its strong dynamics makes the extra quarks 
confined and breaks the accidental PQ symmetry at that 
scale. In this model, a composite NG boson plays the 
role of axion and the PQ mechanism is realized without 
fine-tunings [22], 

Although the above scenario confronts with the cosmo¬ 
logical domain wall problem [23, 24], it can be avoided 
when the PQ symmetry is broken before inflation. In this 
case, however, the axion DM predicts sizable isocurva¬ 


1 See Refs. [19—21] for other mechanisms to explain the origin of 
PQ symmetry. 


ture density perturbations [25-27] due to quantum fluc¬ 
tuations in the axion field during inflation. The resulting 
amplitude of isocurvature perturbations is inconsistent 
with the observation of CMB fluctuations unless the en¬ 
ergy scale of inflation H ln { is smaller than about 10' GeV. 
Such a small energy scale excludes many interesting in¬ 
flation models, such as the Starobinsky model [28] and 
chaotic inflation model [29]. 

In this letter, we investigate a stabilization mechanism 
of the size of extra dimension and propose a scenario 
to suppress axion isocurvature perturbations. 2 We in¬ 
troduce a bulk scalar field that is responsible for the 
stabilization of the size of the extra dimension by the 
Goldberger-Wise (GW) mechanism [38]. Since the bulk 
scalar field is coupled to inflaton on our brane, the size 
of extra dimension changes after inflation. When the size 
of extra dimension during inflation is much smaller than 
the one at present, the PQ symmetry can be drastically 
broken during inflation. As a result, PQ breaking op¬ 
erators give the axion a large mass during inflation and 
suppress undesirable axion isocurvature fluctuations. Af¬ 
ter inflation ends, the size of the extra dimension becomes 
large due to an interaction between inflaton and the bulk 
scalar field, and PQ breaking operators are exponentially 
suppressed. In this scenario, the extra dimension gives 
us mechanisms to explain the origin of the PQ symme¬ 
try and to avoid the axion isocurvature problem. Since 
the Planck scale is proportional to the size of the extra 
dimension, it also changes at the end of inflation. In our 
scenario, dimensionful parameters in the inflaton sector 
should be rescaled by the ratio of effective Planck scales 
during and after inflation. 


II. EXTRA DIMENSION 

Let us consider a 5D spacetime where the 5D manifold 
is R 4 x S/Z 2 just like the RS I model [39] and the metric 


2 See Refs. [30—37] for other mechanisms to avoid the above axion 
isocurvature problem, though they did not mention the origin of 
the PQ symmetry. 
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is written as 

ds 2 = n 2 (f, y)dt 2 - d 2 (f, y)dir 2 - 6 2 (f, y)dy 2 . ( 1 ) 

Here, the coordinate of the extra dimension y extends 
from —1 to 1 and the Z 2 orbifold symmetry is described 
by y —> —y- Two branes are located at the fixed points 
in the S 1 /Z 2 orbifold, i.e., y — 0 and y = 1. 

We introduce matter on branes and a bulk scalar field 
<f> that is responsible for the stabilization of the size of 
the extra dimension by the GW mechanism [38]. The 
action is given by 


( 1 
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where p 5 and g 4 are the determinants of 5D and 4D met¬ 
ric, respectively. The constant n 2 is related to the 5D 
Planck scale M 5 as k 2 = M/ -3 . Hereafter, we rewrite the 
negative bulk cosmological constant A as A = — 6 toq/k 2 . 
£ m ,o and £ m ,i represent the Lagrangians for matter, and 
the energy momentum tensor corresponding to these La¬ 
grangians is described by 


^Latter = ~ (y) diag (p*, ~P*, ~P*, ~P; 0) 

+ ((y - 1)) diag {p, -p, -p, -p, 0). (3) 

6 

Hereafter, we neglect matter on the brane located at y — 
0 and take p* = p* = 0 . 

In this letter, we take a limit of mob -C 1, which implies 
that the extra dimension is almost flat. In this case, the 
solution of the Einstein equation up to the first order of 
p is given as [40, 41] 


n{t,y) = 1 - 
a(t,y) = a(t) 
b(t,y) = b , 
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for 0 < y < 1, where p = up. Throughout this section, we 
consider the case that the radion b stays at its potential 
minimum [see Eq. (14)]. The scale factor a(t) obeys the 
Friedmann equations: 


8t tG 


-P 


-— 47tG (p + p) 

a 


( 7 ) 


( 8 ) 


8ttG~ 


26’ 


( 9 ) 


where a = da/dt. Here, a brane tension on a brane has 
been taken so that the cosmological constant is (almost) 
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TABLE I. Charge assignment for matter fields. 


vanishing in the Friedmann equations. Assuming that 
the VEVs of the bulk scalar field 4> on our brane and the 
other brane are fixed at v± and t7 0 by IA( < I>) and Vo(4 > ), 
respectively, 3 we obtain the static solution of $ as 

$(y) = Ae~ m * by + Be- m * Kl - v) 
v 0 — e- m * b vi 

Y _ g-2m$ b 

^ v\ e~ m ^ b vo 

Y _ g-2m $b ’ 

for 0 < y < 1. Note that <E>, vo, and v\ have mass dimen¬ 
sion 3/2. Substituting this solution into the Lagrangian, 
we obtain the radion potential: 

V r (b) = m* (1 - e~ 2rn * b ) ( A 2 + B 2 ) . (13) 

This implies that the radion stays at 

7?1 

6 ~-In — for v\ > vq, (14) 

770 

and the size of the extra dimension 6 is stabilized. 


( 10 ) 

( 11 ) 

( 12 ) 


III. ORIGIN OF PQ SYMMETRY 

In this section, we briefly review the mechanism to 
explain the origin of the PQ symmetry in the extra di¬ 
mension following Refs. [17, 18]. We introduce a bulk 
hypercolor gauge group SU(N)h and Nf (> 4) pairs of 
fields Q 1 and Q l which transforms under the fundamen¬ 
tal and anti-fundamental representations of SU(N)h , re¬ 
spectively. The first three Q l and Q l of “flavour” indices 
i = 1,2,3 transform under the fundamental and anti¬ 
fundamental representation of SU(3) C , respectively. The 
charge assignment is summarized in Table I. The QCD 
gauge field is assumed to propagate in the bulk while the 
field Q ( Q ) lives on our (the other) brane . 4 In this setup, 
there is an approximate axial symmetry because inter¬ 
actions between the chiral extra-quarks are suppressed 
exponentially. For example, the mass term of 

M q QQ + H.C., (15) 


3 To justify the flat extra dimension, the backreaction of the bulk 

scalar field on the metric should be subdominant. This can be 
achieved when uo,ui k~ 1 [42], 

4 Gauge anomalies on the branes can be cancelled by a Chern- 
Simons term in the bulk [17] . 
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is suppressed as Mq ~ b~ 1 e~ Msb . 

We assume that the hypercolor gauge interaction con¬ 
fines at a scale f a and a chiral condensate develops as 

( 16 ) 

In terms of the SU(N)h gauge theory, the U{Nf)l x 
U(Nf)r flavour symmetry is spontaneously broken to 
the diagonal U(Np)v symmetry and there would be a 
large number of NG bosons in the low energy effective 
theory. However, the flavour symmetry is explicitly bro¬ 
ken by SU( 3) c gauge interactions and SU(3) C charged 
NG bosons acquire masses via SU( 3) c radiative correc¬ 
tions. In addition, there is Z7(l) [SU(N)hI anomaly, so 
that one linear combination of SU( 3) c singlet NG bosons 
obtains a mass through nonperturbative effects. Thus, 
there are (Np — 3) 2 NG bosons in the effective theory 
below the energy scale of f a ■ One of these NG bosons 
that is associated with the following 17(1) symmetry can 
be identified with axion: 

f QfQ*) -t e ia l' i Q\Q i ) for z = 1,2,3 

I QKQ*) -t for i > 4. [ ’ 

This is actually free from U(l) [SU(N)h] 2 anomaly. 
Hereafter we denote this 17(1) symmetry as 17(1)pq. 

The 17(1)pq symmetry has 17(1)pq [S'17(3) c ] anomaly, 
so that the axion acquires a mass through nonperturba¬ 
tive effects after the QCD phase transition [5, 6]. When 
the Hubble parameter becomes comparable to the axion 
mass, the axion starts to oscillate around its minimum [7- 
9] . The present energy density of the axion oscillation is 
calculated as [43, 44] 



where h is the Hubble parameter in units of 
100 km/s/Mpc and 0i n ; is the initial misalignment an¬ 
gle. The domain wall number Apw is equal to N in 
our model because there are N flavours in terms of the 
SU(3) C symmetry. The axion abundance is consistent 
with the observed DM abundance ~ 0.12 when 

the axion decay constant is given as 

fa * 7 x 10 11 GeV x N DW 9r]- 68 , (19) 


where we use Aqcd = 400 MeV. Note that such an 
intermediate scale can be naturally realized in our model 
because the PQ breaking scale f a is determined by the 
strong dynamics of the SU(N)h gauge theory. 

Let us check whether or not PQ breaking operators 
are sufficiently suppressed to explain the smallness of the 
strong CP phase. The operator of Eq. (15) gives the 
axion a mass of order \jMQf a at a minimum that is 
generally different from the one determined by the origi¬ 
nal vacuum angle 9q. Thus, it may induce a shift in the 
strong CP phase as 


Ad ~ 


( VWa \ 
y rn a J 


( 20 ) 


where m a is the axion mass at the low energy: 

m u m d m^f n ^ ^ 

Here, f n is the pion decay constant, and m u , mj, and m 
are the masses of up quark, down quark, and pion, respec¬ 
tively. Requiring A 9 < 0 o bs ~ 1CU 10 and using Mq ~ 
b~ 1 e~ M5b , we find that M 5 b > 150 + 31n(/ Q /10 12 GeV) 
is sufficient to solve the strong CP problem. 

However, operators involving either Q or Q are ex¬ 
pected to be suppressed only by powers of M^ 1 on each 
brane. In the case of N = 3, for example, we can write 
the following terms on the branes: 

J d*xV=r<$s {QQQ? 

+ j^xV^I^(QQQ) 2 + n. c ., ( 22 ) 

where i/q and uq are coupling constants. These terms 
also give a mass to the axion and may shift the strong 
CP phase as 

A9-5X10- 1 3 W 5 (_4_) 11 (A.) ‘V) 

Noting M 5 = Mpi/V 2 &M 5 < Mpi/18, this is marginally 
consistent with the present upper bound on the strong 
CP phase for yQUQ < 10 -3 . Although the bulk scalar 
field $ and radion b can interact with Q and Q on 
each brane, these contributions are also suppressed in the 
same way. Note that A 9 is much more suppressed when 
we consider N > 5. In the case of IV = 5 and Np > 5, 
for example, the lowest dimension operator is written as 
{QQQQQ) 2 /Ml 1 and the shift in the strong CP phase is 
at most of order 10 _5 ° for f a = 3 x 10 12 GeV. Note that 
we have assumed that light colored particles other than 
Q and Q are absent in the bulk and on the branes. This 
ensures that terms like Eq. (22) are the lowest dimension 
operators, so that explicit breakings of the PQ symmetry 
is sufficiently suppressed to solve the strong CP problem. 


IV. SUPPRESSING ISOCURVATURE 
FLUCTUATIONS 

When the PQ symmetry is broken during inflation 
and the mass of axion is much less than the Hubble 
parameter, the axion acquires sizable quantum fluctua¬ 
tions and predicts isocurvature perturbations which can 
be observed in CMB fluctuations [25-27]. However, the 
observation of CMB fluctuations reveals that they are 
predominantly adiabatic and puts an upper bound on 
the amount of isocurvature perturbations. This upper 
bound implies that the energy scale of inflation has to be 
smaller than of order 10 7 GeV [45]. Such a small energy 
scale severely restricts inflation models. This is called 
the axion isocurvature problem. 
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In our model, the mass of axion is suppressed exponen¬ 
tially as e~ M ^ b (see Eq. (15)). This implies that if the 
size of the extra dimension b during inflation is smaller 
than its present value by more than one order of magni¬ 
tude, the mass of axion can be larger than the Hubble 
parameter during inflation and axion fluctuations can be 
suppressed. This scenario can be realized in the following 
way. Let us consider a scalar field ip (such as an inflaton 
or a waterfall field) which lives on our brane and whose 
VEV changes from 0 to ipo after inflation. We assume 
that p and the bulk scalar field $ have an interaction 
term such as 


Vi($) 


A ' + Ai, (24) 


4 Mf 


where Ai and A; n t are dimensionless parameters. Then 
the VEV of $ on our brane changes from iq to v[ = 
i Jv\ + (Ai n t/ 2 Ai)M 5 (^Q after inflation, which means that 
the size of the extra dimension changes after inflation 
via Eq. (14). As a result, the mass of axion can be much 
larger than its present value and can be larger than the 
Hubble parameter during inflation. For example, when 
M 5 b' n{ ~ 10 during inflation, the axion mass is as large 
as m a ~ y/M Q f a ~ 10 12 GeV. 

In order to solve the strong CP problem, the axion 
mass has to be suppressed and the size of the extra di¬ 
mension should satisfy M 5 6 now > 150+3 ln(/ o /10 12 GeV) 
at the present epoch. In the case of v\/v 2 = 1.1, 
v'i/v 2 = 5, and m$/M 5 = 0.01, for instance, we ob¬ 
tain desired values: M$b' ni ~ 9.5 during inflation and 
M 5 b now ~ 161 at present. Note that the axion VEV 
during inflation is generally different from the one at the 
present epoch. This is because the former one is deter¬ 
mined by the phase of the mass term of Eq. (15) while the 
latter one is determined by the original vacuum angle 0 O ■ 
Therefore the axion starts to oscillate around the latter 
minimum at the QCD phase transition and its abundance 
is given by Eq. (18). 


V. PREDICTIONS 
A. inflation 

In our scenario, the extra dimension is almost flat and 
its size is of order the 5D Planck length so that high 
scale inflation can be realized. Although such a small ex¬ 
tra dimension cannot be directly accessible at collider ex¬ 
periments, predictions of CMB fluctuations may be mod¬ 
ified in our scenario. This is because the Planck scale dur¬ 
ing inflation Mp" f is different from the one in the present 
epoch Mp° w as Mp“ f = <V/ 2 Mp° w (see Eq. (9)), where 
S = b int /b n ™ (= O( 10 -1 )). 

For simplicity, we consider the case that the size of 
the extra dimension is constant during the time that the 
CMB modes exit the horizon. Since curvature pertur¬ 
bations are conserved outside the horizon, the result¬ 
ing CMB spectrum is easily calculated by rescaling the 


Planck scale. During inflation, the Hubble parameter is 
given by [46] 


H 2 



1 

3(5Af 2 ! 


V. 


The spectral index is calculated from n s = 1 
where the slow roll parameters are given by 

e =^(f) 2 

V = SMpY (ypj , 


(25) 
6 e + 2r /, 


(26) 

(27) 


where the prime denotes a derivative with respect to in¬ 
flaton. These results imply that slow roll inflation re¬ 
quires a steeper potential than ordinary inflation mod¬ 
els. The amplitude of scalar perturbations of the metric 
is given by 


Ak ~ <5 2 M^24tt 2 ’ 

and the tensor-to-scalar ratio is given by r = 16e. It 
is clear that once all dimensionful parameters in the 
inflaton sector are rescaled by a factor of <5 1//2 (e.g., 
+inf —> (5 1 ^ 2 +inf and V —> 5 2 V), we reproduce standard 
predictions of inflation. This is because the rescale of the 
Planck scale can be performed by a conformal transfor¬ 
mation, which does not change the values of the dimen¬ 
sionless parameters e, 77 , and A#.. 

We assume that the VEV of ip in Eq. (24) changes from 
0 to ipo after inflation. This is naturally realized when 
we consider a hybrid inflation model and identify ip as 
a waterfall field. Denoting inflaton as J, we write the 
potential in the inflaton sector as 

V in{ =p 2 I 2 + pip A + Vi(I). (29) 

The inflaton / has a large VEV during inflation 
and slowly rolls toward the critical point I c = 
(Aint+i/KM 5) 1 / 2 due to an inflaton potential V/(7). 
When the VEV of the inflaton decreases to J c , the water¬ 
fall field ip starts to oscillate around the true minimum 
and the inflation ends. At the same time, the VEV of 
the bulk scalar field $ on our brane changes from v\ to 
v[ and the size of the extra dimension changes from b lnt 
to 6 now . 

Finally, we should note that the oscillations of $ and 
the radion b may be induced after inflation due to the 
change of their potential minimum. When we introduce 
an interaction of QNrNr on our brane, where TV# is 
the right-handed neutrino, $ can decay into right-handed 
neutrinos and then the right-handed neutrinos decay into 
standard-model particles. Since the radion b couples with 
SM particles via Planck-suppressed operators, it can also 
decay into radiation. In our model, the mass of the radion 
is given by 


ml = 


3Mpj 


(& now ) 2 m|i 


W - vz 


(30) 
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after its kinetic term is canonically normalized. This is 
of order 10 13 GeV for Vi,v 2 = k _1 /100, so that the re¬ 
heating temperature can be as large as 10 10 GeV even if 
they dominate the Universe. Such a high reheating tem¬ 
perature is consistent with the realization of the thermal 
leptogenesis [47]. 


B. dark radiation 


In our model, there are (Np — 3) 2 — 1 massless NG 
bosons in addition to the axion, so that they may con¬ 
tribute to the energy density of the Universe as dark ra¬ 
diation [48-50] . Since they interact with the SM particles 
via interactions suppressed by some powers of f a , their 
decoupling temperature is roughly an order of magni¬ 
tude below f a (see e.g., Refs. [51, 52]). Therefore if the 
maximal temperature after inflation is lower than f a and 
the reheating temperature can be as large as of order 
/ a /10, the NG bosons may be thermalized after infla¬ 
tion without restoring the PQ symmetry. In fact, such 
a high reheating temperature is favoured to realize the 
thermal leptogenesis [47]. Once the NG bosons and ax¬ 
ion are thermalized and decoupled at a high temperature, 
they contribute to the energy density of the Universe as 
dark radiation. Its amount is conventionally expressed by 
the effective neutrino number and the result is calculated 
as [52] 


N eS ~ N^ m) + 0.027 x (N f - 3) 2 , (31) 


where (~ 3.046) is the SM prediction. The present 

constraint on N e g is N eS = 2.99 ± 0.39 (95% C.L.) [53, 
54], which implies that Np should be smaller than or 
equal to 6. The ground-based Stage-IV CMB polariza¬ 
tion experiment CMB-S4 will measure the effective neu¬ 
trino number with a precision of AiV e ff = 0.0156 within 
one sigma level [55] (see also Ref. [56]). The number of 
flavours can be measured indirectly via the observation 
of the amount of dark radiation. 5 


VI. DISCUSSION AND CONCLUSIONS 

We have proposed a QCD axion model in a 5D space- 
time to explain the origin of the PQ symmetry and 
avoid the axion isocurvature problem. Each extra quark- 
antiquark pair lives on branes separately so that they 
have an accidental PQ symmetry [17, 18]. We assume 
an interaction between inflaton and a bulk scalar field 
that is responsible for the stabilization of the size of the 
extra dimension. As a result, the size of the extra di¬ 
mension during inflation can be different from the one at 
present and PQ breaking operators can be efficient dur¬ 
ing inflation. We have shown that the axion mass from 
PQ breaking operators can be larger than the Hubble 
parameter during inflation to suppress axion isocurva¬ 
ture fluctuations while the PQ symmetry is sufficiently 
preserved after inflation to solve the strong CP problem. 
Therefore, the extra dimension explains not only the ori¬ 
gin of the PQ symmetry but also suppression of axion 
isocurvature fluctuations. 

Since the size of the extra dimension changes after in¬ 
flation, all dimensionful parameters in the inflaton sector 
should be rescaled by a certain factor compared with the 
ones in ordinary inflation models. In addition, our model 
may predict a sizable amount of dark radiation which 
would be detected in the near future. Finally, we com¬ 
ment on our assumption that the bulk scalar $ and the 
ration b stay at their minimum. This is actually justified 
during and well after inflation because of the redshift ef¬ 
fect. The dynamics of $ and b just after inflation might 
be complicated and will be discussed elsewhere. 
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Although the amount of dark radiation for Nf = 6 is exactly 
the same with the one predicted in Ref. [52] , we could distinguish 


between these models via the observation of an exotic kaon decay 
process. 





6 


[1] C. A. Baker, D. D. Doyle, P. Geltenbort, K. Green, 
M. G. D. van der Grinten, P. G. Harris, P. Iaydjiev and 
S. N. Ivanov et al., Phys. Rev. Lett. 97, 131801 (2006) 
[hep-ex/0602020]. 

[2] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 
(1977). 

[3] R. D. Peccei and H. R. Quinn, Phys. Rev. D 16, 1791 
(1977). 

[4] S. Weinberg, Phys. Rev. Lett. 40, 223 (1978); F. Wilczek, 
Phys. Rev. Lett. 40, 279 (1978). 

[5] G. ’t Hooft, Phys. Rev. Lett. 37, 8 (1976). 

[6] G. ’t Hooft, Phys. Rev. D 14, 3432 (1976) [Erratum-ibid. 
D 18, 2199 (1978)]. 

[7] J. Preskill, M. B. Wise and F. Wilczek, Phys. Lett. B 
120, 127 (1983). 

[8] L. F. Abbott and P. Sikivie, Phys. Lett. B 120, 133 
(1983). 

[9] M. Dine and W. Fischler, Phys. Lett. B 120, 137 (1983). 

[10] S. B. Giddings and A. Strominger, Nucl. Phys. B 307, 
854 (1988). 

[11] S. R. Coleman, Nucl. Phys. B 310, 643 (1988). 

[12] G. Gilbert, Nucl. Phys. B 328, 159 (1989). 

[13] T. Banks and N. Seiberg, Phys. Rev. D 83, 084019 (2011) 
[arXiv:1011.5120 [hep-th]]. 

[14] M. Kamionkowski and J. March-Russell, Phys. Lett. B 
282, 137 (1992) [hep-th/9202003]. 

[15] S. Ghigna, M. Lusignoli and M. Roncadelli, Phys. Lett. 
B 283, 278 (1992). 

[16] B. A. Dobrescu, Phys. Rev. D 55, 5826 (1997) [hep- 
ph/9609221], 

[17] K. I. Izawa, T. Watari and T. Yanagida, Phys. Lett. B 
534, 93 (2002) [hep-ph/0202171]. 

[18] K. I. Izawa, T. Watari and T. Yanagida, Phys. Lett. B 
589, 141 (2004) [hep-ph/0403090]. 

[19] R. Holman, S. D. H. Hsu, T. W. Kephart, E. W. Kolb, 
R. Watkins and L. M. Widrow, Phys. Lett. B 282, 132 
(1992) [hep-ph/9203206]. 

[20] S. M. Barr and D. Seckel, Phys. Rev. D 64, 123513 (2001) 
[hep-ph/0106239]. 

[21] K. Harigaya, M. Ibe, K. Schmitz and T. T. Yanagida, 
arXiv:1505.07388 [hep-ph]. 

[22] J. E. Kim, Phys. Rev. Lett. 43, 103 (1979); M. A. Shif- 
man, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 
166, 493 (1980). 

[23] Y. B. Zeldovich, I. Y. Kobzarev and L. B. Okun, Zh. 
Eksp. Teor. Fiz. 67, 3 (1974) [Sov. Phys. JETP 40, 1 
(1974)]. 

[24] P. Sikivie [ADMX Collaboration], Phys. Rev. Lett. 48, 
1156 (1982). 

[25] M. Axenides, R. H. Brandenberger and M. S. Turner, 
Phys. Lett. B 126, 178 (1983). 

[26] D. Seckel and M. S. Turner, Phys. Rev. D 32, 3178 
(1985). 

[27] M. S. Turner and F. Wilczek, Phys. Rev. Lett. 66, 5 
(1991). 

[28] A. A. Starobinsky, Phys. Lett. B 91, 99 (1980). 


[29] A. D. Linde, Phys. Lett. B 129, 177 (1983). 

[30] A. D. Linde and D. H. Lyth, Phys. Lett. B 246, 353 
(1990). 

[31] A. D. Linde, Phys. Lett. B 259, 38 (1991). 

[32] S. Kasuya, M. Kawasaki and T. Yanagida, Phys. Lett. B 
409, 94 (1997) [hep-ph/9608405]. 

[33] S. Kasuya, M. Kawasaki and T. Yanagida, Phys. Lett. B 
415, 117 (1997) [hep-ph/9709202]. 

[34] S. Folkerts, C. Germani and J. Redondo, Phys. Lett. B 
728, 532 (2014) [arXiv:1304.7270 [hep-ph]]. 

[35] K. S. Jeong and F. Takahashi, Phys. Lett. B 727, 448 
(2013) [arXiv: 1304.8131 [hep-ph]]. 

[36] M. Kawasaki, T. T. Yanagida and K. Yoshino, JCAP 
1311, 030 (2013) [arXiv: 1305.5338 [hep-ph]]. 

[37] K. Nakayama and M. Takimoto, arXiv:1505.02119 [hep- 
ph[. 

[38] W. D. Goldberger and M. B. Wise, Phys. Rev. Lett. 83, 
4922 (1999) [hep-ph/9907447]. 

[39] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 
(1999) [hep-ph/9905221]; Phys. Rev. Lett. 83, 4690 
(1999) [hep-th/9906064]. 

[40] C. Csaki, M. Graesser, L. Randall and J. Terning, Phys. 
Rev. D 62, 045015 (2000) [hep-ph/9911406]. 

[41] J. M. Cline and H. Firouzjahi, Phys. Lett. B 495, 271 
(2000) [hep-th/0008185]. 

[42] O. DeWolfe, D. Z. Freedman, S. S. Gubser and A. Karch, 
Phys. Rev. D 62, 046008 (2000) [hep-th/9909134], 

[43] O. Wantz and E. P. S. Shellard, Phys. Rev. D 82, 123508 
(2010) [arXiv:0910.1066 [astro-ph.CO]]. 

[44] M. Kawasaki, K. Saikawa and T. Sekiguchi, Phys. Rev. 
D 91, no. 6, 065014 (2015) [arXiv:1412.0789 [hep-ph]]. 

[45] P. A. R. Ade et al. [Planck Collaboration], 
arXiv: 1502.02114 [astro-ph.CO]. 

[46] D. I. Kaiser, Phys. Rev. D 52, 4295 (1995) [astro- 
ph/9408044], 

[47] M. Fukugita and T. Yanagida, Phys. Lett. B 174, 45 
(1986). 

[48] M. S. Turner, Phys. Rev. Lett. 59, 2489 (1987) [Erratum- 
ibid. 60, 1101 (1988)]. 

[49] K. Nakayama, F. Takahashi and T. T. Yanagida, Phys. 
Lett. B 697, 275 (2011) [arXiv: 1010.5693 [hep-ph]]. 

[50] S. Weinberg, Phys. Rev. Lett. 110, no. 24, 241301 (2013) 
[arXiv:1305.1971 [astro-ph.CO]]. 

[51] E. Masso, F. Rota and G. Zsembinszki, Phys. Rev. D 66, 
023004 (2002) [hep-ph/0203221], 

[52] M. Kawasaki, M. Yamada and T. T. Yanagida, 
arXiv: 1504.04126 [hep-ph], 

[53] E. Aver, K. A. Olive, R. L. Porter and E. D. Skillman, 
JCAP 1311, 017 (2013) [arXiv: 1309.0047 [astro-ph.CO]]. 

[54] P. A. R. Ade et al. [Planck Collaboration], 
arXiv: 1502.01589 [astro-ph.CO]. 

[55] W. L. K. Wu, J. Errard, C. Dvorkin, C. L. Kuo, 
A. T. Lee, P. McDonald, A. Slosar and O. Zahn, Astro- 
phys. J. 788, 138 (2014) [arXiv: 1402.4108 [astro-ph.CO]]. 

[56] K. N. Abazajian et al. [Topical Conveners: K.N. Abaza- 
jian, J.E. Carlstrom, A.T. Lee Collaboration], Astropart. 
Phys. 63, 66 (2015) [arXiv:1309.5383 [astro-ph.CO]]. 


